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Abstract 

We propose a new method to describe the interacting bose gas at zero temper- 
ature. For three-dimensional system the correction to the ground- state energy 
in density is reproduced. For two-dimensional dilute bose gas the ground- state 
energy in the leading order in the parameter | Ino^pl^^ where a is a scattering 
length is obtained. 
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At present time two-dimensional models attract much attention in connection with 
the problems of high- Tc superconductivity and fractional quantum Hall effect. In this 
context the description of the two-dimensional system of bosons may be important. For 
instance the system of planar fermions in the magnetic field is equivalent to the system 
of the interacting bosons with the additional long-range Chern-Simons interaction. 
The problem of the hard core bosons on a lattice at high density is closely related to 
the description of different strongly correlated electronic systems. The conventional 
methods of the description of the interacting bose gas are inapplicable in both of the 
above mentioned problems. In this context the new ways to treat the system in the 
two different physical limits where the perturbation theory is possible are of interest. 

Various methods to describe the system of bosons with the pairwise interaction 
were introduced which in one way or another consisted in summing up an infinite 
subset of the terms in the perturbation series. The example of this procedure is the 
Bogoliubov's method which is correct in the high density limit which means that 
the range of the potential is much larger than the average particle spacing. In the 
opposite limit of the dilute bose gas the perturbation theory does not work which 
manifests in the divergence in the expression for the energy for the (5-function type 
potential. Various modifications of this procedure using the pseudopotential or the 
diagrammatic expansion methods, for example, were proposed . However the reason 
why the procedure of ref . ^ can be applied when the two-body potential is not small and 
the wave function is not close to the unperturbed one is obscured. The generalization 
of this methods to the 2D system is not straightforward due to the behavior of the 
modified scattering amplitude which approaches zero in the low energy limit. In fact 
the methods can be considered as an application of the Bogoliubov's approach to a 
system with the parameters chosen in such a way that both this approximation is valid 
and at the same time the density p is small in a sense ap^^^ « 1 (a is the scattering 
length). Although this can be a basis for the solution of 2D problem the estimate of 
the accuracy of the approximation [|I| for a given potential is required. 
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In the present paper we suggest a new method to describe the system of bosons for 
the two hmiting cases and investigate the ground state properties of 2D bose gas at zero 
temperature. Our method is closely related to the approach proposed by Lieb which 
makes use of the distribution functions related to the ground state wave function. We 
use the decomposition of the logarithm of the wave function into a series over the n- 
particle functions. We argue that at small density the proposed expansion corresponds 
to the expansion of energy in the small parameter. Although in practice it is difficult to 
solve the corresponding system of equations beyond the leading order approximation 
which in the framework of our approach corresponds to the wave function of the Jastrow 
form, the convincing arguments in favor of the validity of this procedure in the leading 
order can be found. Of course the distribution functions can be calculated using the 
wave function of the Jastrow form. In this way the connection with the approach of 
ref.0 can be made. 

The equation for the energy Eq and the wave function $(xi, . . . xtv) of the ground 
state for a system of particles in the volume V interacting with the two-body po- 
tential U{x) has the form 

{-Y.d} + 2Y,U{x,M^ = Eo^, z,j = l,...Ar, (1) 

where di = d/dxi, and we denote by Xj the d- dimensional space vector 

of i-th particle throughout the paper. The ground state wave function ^{xi, ...xn) is 
a symmetric and positive function of its arguments. According to (see also 0) one 
can seek for the function $ in the form 

^{xi, . . . xj\[) = exp S{xi, . . .xjsf). 

We observe that the function S{xi, . . .xn) can be expanded as follows 

S{xi,...xn) = J2^2(^ii) + Sai^i, Xk, xi) + . . . , (2) 

i<j i<k<l 
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where SJ the symmetric functions depending on the relative interparticle 

distances and subjected to the constraints 

J dXnSniXi, . . .Xn) = 0, U > 2. (3) 

Integrating Eq.(0) over the coordinates of N — n particles subsequently for n = 2,3 . . . 
and using the condition (j^) one can show that the expansion (0) is an irreducible one 
which means that for a given function S{xi, . . . xn) the set of the functions Sn is unique. 
Substituting Eq.(||) into the equation (|l]) we obtain the equation 

Eo = Y. ^(^u) - E dS{x,k)dS{xa) + ..., (4) 

F{x) = -d^S{x) - dS{x)dS{x) + U{x), 

where the terms depending on the pair function 5*2 (x) = S{x) only are indicated 
explicitly. Integrating Eq.(||) over the coordinates n + 1, . . . and using Eq.(|]) we 
obtain the n-particle equation for the functions Sn- In general the system of the 
equations cannot be solved and the assumptions on the higher order Sn - functions 
are required. For the dilute bose gas the perturbation theory in the small parameter 
can be used. Although the calculations in the high orders are too complicated one can 
argue that the expansion of the energy in the small parameter corresponds just to the 
expansion (^. Here we will demonstrate it in the lowest order which corresponds to 
the pairwise wave function. It is plausible that the same is valid for the higher order 
terms. Suppose for a while that only the function S{x) is not equal to zero. Integrating 
the equation (^) over the coordinates 1, ... we obtain for the parameter a defined as 

Eo = Npa 

(p = N/V, Eq/V = ap^) the equation 

A^ - 1 f 

a = — — — J dxF{x) (5) 
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which is equivalent to a = / dx{U — {dS)"^) in the infinite volume limit. Integrating 
Eq.@ over the coordinates 3, . . . and using Eq.(|) we obtain the equation 

d'^S{x) + dSdS{x) - U{x) - pdl j dx^S{xi^)S{x2z) = ~a/V. (6) 

where the notation x = xu is used and the terms of order ~ 1/V are retained. The 
right- hand side is equal to zero in the limit V ^ oo. We assume the periodic boundary 
conditions so that the integral over the total derivative is zero. Due to the condition (^) 
the contribution of S3 to Eq.(|) is of the form p J dxidx2{diS{123)y and the example 
of its contribution to Eq.(||) is p / dx^d^S (123) d^S (13). One can use the solution of 
Eq.(l^) to estimate the function S3 with the help of the three-particle equation and 
show that these integrals are suppressed. We will show below that it is the smallness 
of S{x) at the distances ~ p^^^^ (d = 2, 3) that allow one to neglect the 5*3 function in 
the lowest order approximation. The equation reproduce the correction in density 
[0] to the ground- state energy for 3D system and predicts the leading order term for 
2D system. 

Note also that Eq.(|^) can be regarded as an equation for the trial variational wave 
function of the Jastrow form. In this case it follows from the form of the solution of 
Eq. (1^) that the variational energy is given by Eq. with the accuracy up to the higher 
order terms in the expansion parameter. 

If the second term in the left-hand side of the equation (|^) can be neglected then 
the equation can be easily solved using the Fourier transformation: 

pSl - k^Sk -Uk = 0, Sk = J dxe'^'^Six). 

Substituting the solution of this equation Sk to the expression for the energy 

a = f/o + j dx{dS{x)f = Uo + k^Sl, 

{Uq = J dxU{x)) we get 

a = f/o + ^ ((A:' + mpk'') - fc' - 2Uu(^ , (7) 
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where the notation = J d'^k /{2ttY is used. Eq.(^ is the Bogohubov's expression for 



the energy. One can estimate the accuracy of Eq.(^ for a given two-body potential 
U{x). The corrections are determined by the function S{x) which should be small in 
order the approximation (|^) to be valid. For instance in 2D for the potential which 
is a repulsive square well in momentum space, Uk = u, k < b^^, Uk = 0, k > b~^, 
the energy is: (i) at upb'^ » 1, a = u{l - 4^^), 5'(0) = -{u/ pb'^Y^'^ /2tc; (ii) at 
upb'^ « 1, a = u{l - ^\n{upb^e^/^)), S{0) = -l/^npb'^. Briefly, the potential should 
be relatively shallow compared to its width and the spatial range should be much larger 
than the average particle spacing. 

Let us proceed with the solution of the equation in the low density limit. The 
particle density is assumed to be small so that the dimensionless parameter ap^^'^ « 
1. In this formula a is the scattering length for the potential U and p-^/'^ is the 
average interparticle distance in d- dimensional space. The expansion parameters are 
respectively |lnpa^|~^ and o?l'^p^l'^ in two and three dimensions. 

1. Three-dimensional system. 

Let us begin with the case of three spatial dimensions d = 3. Define Sq{x) as the 
solution of the equation (||) at p = 0: 



Here a is the unknown parameter and the periodic boundary conditions are implied. 
In terms of the function (^{^[x) = exp(5'o(x)) keeping the volume finite and taking into 
account the periodic boundary conditions we find that 



is proportional to the scattering length for the potential U{x). For the hard-sphere 
potential a is the radius of the potential. In the infinite volume limit the equation 
has the form {—d^ + ?7)0o(x) = which is to be supplemented by the boundary 



d^So + dSodSo -U = -a/V. 
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a = 




6 



condition (j)o{r) 1 at r oo {r = \x\). At the distances much larger than the 
range of the potential the solution is (f)o{r) = 1 — a/r. For the energy we have a = 
a + J dx^S&^S — SqO'^Sq) or, equivalent ly, 

a = a- UeSl-k^Sl). (9) 

Jk 

Make use of the Fourier transformation in the equations and (^). In the region 
r » a where the condition Sq{x) << 1 is satisfied we have 5*0 (r) = —a/r and the 
main contribution to the integral / dxldSo)"^ comes from the integration over the region 
of the small r ~ a. Up to the corrections in p the same is true for the analogous integral 
for S{x). 

Therefore the Fourier transform of the function U — [dS^, and the same function 
for So, are independent of the momentum k a.t k « 1/a. At these values of k, Sok = 
—a/k"^ and Sk is given by the solution of the equation pk^Sl—k^Sk—a = 0. Substituting 
these functions into Eq.(P) we obtain 



(10) 



The integral in this expression converges at large k and is saturated at the values 
k ~ {apY^"^ << 1/a which justifies our assumption that / dxe'^^^{U — (dSY) = a. The 
corrections to this formula does not change the result in the approximation considered. 
Evaluating the integral in Eq. (|l^) we obtain the well known result for the expansion 
of the ground state energy in powers of density in 3D: 

2. Two-dimensional system. 

Let us proceed with the investigation of two-dimensional system. In 2D it is easy 
to find the solution of Eq. (P) and determine the energy a as the function of the density 
p with the accuracy up to the terms of the higher order in the small parameter a ~ 
I In pa'^\~^. 



As in three-dimensional case the solution of Eq.(^) can be represented in the form 



where 5"^ is the Fourier transform of the function S{x) and the function cxfc is defined 
by 

ak = J dxe'''^ {U{x) - dSdS{x)) . 

At the sufficiently small k we have (Xfc = a while at the larger k, ak could depend on the 
momentum k. However it can be shown that in fact with the accuracy of order O(a^) 
one can take ak — a for all << a~^. We are interested in the behavior of the function 
S{x) in the region a « r « (ap)^^/^. As in 3D (ap)^^/^ is the characteristic length, 
the kind of correlation length in the problem. In this region S{x) can be calculated as 

SW--i/e-ffl + ^V''-l) (12) 



2pJk \\ P 
Evaluating this expression we find that 

S{x) = -a/An + (a/27r) In ((ap)^/V) (13) 

which is valid in the indicated region with the accuracy up to the terms of order of 
a^. One can see from Eg . (|T^) that the function S{x) decreases quickly at the distances 
much larger than the correlation length (ap)^^/^: 

Now using the equation (|T3|) one can estimate the function ak- 

- a\ < — In —= + 0[a ). 
In yUip 

This relation justifies the initial assumption about the behavior of ak- The other way 
to obtain the solution (|13D is to represent the equation (||) in the form 



{d^ - f/(x))0(x) = p<j){x)dl J dxsS{x,3)S{x23) - {a/V)<f){x) (14) 



where = expS'(x). The Fourier transform of the function U{x)(f){x) does not 

depend on the momentum at k << 1/a. In the large distance region where S{x) << 1 
we can expand the function (f){x) ~ 1 + S{x) and neglect the terms of order S'^(x). Then 
the equation ( p!4D has the solution of the same form as Eq.(|13[). The only difference is 
that the parameter a is replaced by the value of the integral / dxU{x)(t>{x). To calculate 
it one has to integrate both sides of the equation (|1^) over the space vector in the finite 
volume V: 

dxU {x)(j){x) = a — p I dx{l — (j){x))xix) (15) 



In general the limit ^ — > oo of integral of some function over the volume does not 
coincide with the integral of the limiting function, because in the finite volume the 
function may have the asymptotic const. /V. This is indeed the case for / dxx{x) 
which is zero in the finite volume. However one can substitute the product of the 
limiting functions (1 — 0(x))x(x) into Eq.(|T5D since both functions approaches zero at 
r oo. Alternatively, one can use the equation @ directly in the infinite volume 
limit. The function x(r) can be evaluated using Eq.(|12|): 



1 /a^/^\ 1 

X(r) - aS{r), {a « r « (ap)-^/^); x{r) ^ — [yj^j » M"^^^) 



The integral in the right hand side of the equation (|T5D is determined by the long- 
distance region {S{x) « 1) and can be estimated as O(a^). Thus the solution of 
Eq.d) function of the parameter a is found. 

On the other hand at the sufficiently small r the function (f){x) can be found ap- 
proximately as the solution of the equation (—9^ -|- U{x))(f){x) = due to the fact that 
the right hand side of Eq.(p!4|) is ~ p. The solution is (poir) = Cln(r/a) where C is 
an arbitrary constant and a is the scattering length for the potential U (r) (the region 
r > > a is implied) . In two dimensions the scattering length is defined by the behavior 
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of the scattering amplitude at low energy which (for our equation) is given by 

^"(^) = 1 . lol -1^ + ^(^")' 

m Ka/z + 7 — iTi / 2 

where \ f{k)\ = (27rfc)^/^|/(fc) | is the modified scattering amplitude, k is the momentum 
and 7 = 0.1159 is the Euler constant. For the 2D hard-sphere potential the scattering 
length a is equal to the radius of the potential. The correction 6(f){r) ~ apr^ to the 
solution of the homogeneous equation can be easily estimated using Eq.(^). The 
correction is ~ 1 at the distances of order of the correlation length (ap)~^/^ however 
it is small (50(r) ~ a) at the distances r ~ p^^l'^ . At r ~ p^^/^ the equation (|13D is 
still valid. Comparing the solution given by Eq.(|T^) with the function 0o(r) we get the 
relation 

Cln(r/a) = 1 - a/47r + (a/27r) In ((ap)-^/V) (16) 

which should be valid in the region a « r « (ap)^^/^ with the logarithmic accuracy. 
In particular at r ~ p^^/^ it is valid with the accuracy up to the terms of order ~ a. 
Thus in the leading order in the small parameter | Inpa^l"^, C = a/2T[ and the energy 
is 

-^.-0{^V ,17) 



I In pa^l \y I Inpa 
Eg . ([T7|) is our final result for the two- dimensional system. 

Naively, from the three-particle equation the function 5*3 can be estimated by an 
order of magnitude as 5'(1,2,3) ~ S'(13)S'(23) + permutations. Then the contribu- 
tion of ^3 to Eq.(^) at the distances ~ p~^/^ which is determined by the integral 
p / dx^S{l?>){didz — 9|)S'(123) (see Eq.(D) is suppressed due to the smallness of the 
function S{x) at these distances. Actually the three-particle equation can be solved in 
the leading order in the expansion parameter (that means that the terms ~ >S'2S'3 in this 
equation should be neglected) in the momentum representation and the corresponding 
integral in the momentum space can be estimated (see Appendix). 

At the same time at short distances where the two particles interact strongly the pair 
function is not expected to be very different from the solution of the two-body problem 
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whether or not the function S3 is taken into account. In the other words although the 
function S3 is not necessary small at short distances our estimate of the right hand 
side of Eq.([T^) as well as the equation (|TB|) are valid by an order of magnitude. Note 
that the "momentum" corresponding to the two-particle equation is of order of the 
correlation length which characterizes the screening of the pair wave function due to 
the other particles. Hence the result ( p!7D has a simple physical interpretation. Namely 
the energy of the two particles interacting with the potential U located in the two- 
dimensional volume V is ^ f /V where the modified scattering amplitude / should be 
normalized at /c ~ p^/^. Multiplying the result by the number of pairs we obtain the 
result ([T7|) . In three dimensions that corresponds to the first term in Eq . (|TT|) . 

The other approach to the description of the dilute bose gas based on the equations 
for the distribution functions related to the function $ was proposed in ref. 0]. We will 
show below that the assumptions made in ref. ||^ can be justified in the framework of 
our approach. Let us describe briefly this method. The n-particle distribution function 
is defined by 

gn{xi, ...Xn) = ^"V" J dxn+i ■ ■ ■ dxN^{xi, ...Xn) (18) 
Z = y dxi . . . dxN^{xi, . . . Xn)- 
The energy is related to the pair distribution function g2{xi2) = g{x) by means of 

- 1 r 

a = J dxU{x)g{x). (19) 

Using Eq.(|ll) we obtain the equation for g{x): 

(-d^ + U{x)) g{x) = \Eog{x) -p [ dx3g3il, 2, 3)[/(x23) (20) 



2 

1 (A^-2)(Ar-3) 



dxsdx4g4^{l, 2, 3, 4)f/(a;34). 



2 

The first and the third terms in the right- hand side are of order ~ while the left- 
hand side is of order ~ 1. In the leading order the cancellation of these terms take 
place and the corrections ~ 1/iV are important. The approximation for the distribution 
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functions and used in ref. to evaluate the right hand side of Eq.(pOD is the 
superposition approximation: 

n 

gn{xi,...Xn) = Y[gixij) (21) 

i<j 

From physical considerations it is clear that in the limit V oo the function g{x) — >■ 
g{x) while for the finite volume it can be found using the set of the equations 

(Xi, ...Xn) = Vgn-l{Xl, . . 

For instance, for the four- particle function g^, substituting the ansatz ( pID into the 
formula / dx^dx^ g^i^l, 2, 3, 4) = V'^g{xi2) we find that the corresponding pair function 
is 

~g^'\x) = g{x) (l - ^ / rfX3/(xi3)/(x23)) , 

where the function f{x) = 1 — g{x) can be taken in the infinite volume limit. These 
corrections are not important in the second term of the right- hand side of Eq. (pOD . 
Substituting the expressions for the functions gs and g4 into the equation (pO|) and 
taking the limit — oo we get an equation for the pair distribution function g{x). 
Assuming that the function /(x) is small and retaining the terms of the leading order 
in /(x) we get 

{-d^ + U{x)) g{x) = 2apf{x) - ap^ J dx^f {x is) f{x 23) (22) 

This equation is valid for those r where the function f{x) is small (/(x) << 1). 

The equation ( ^2]) can be solved in the same way as Eq.(|^). For instance the 
solution of Eq. (^2]) for 2D problem in the region where /(x) << 1 is 

g{x) = 1 + a/ATT + (a/27r) In ((ap)^/V) (23) 

for r « (ap)~^/^. Repeating the arguments leading to Eq.(|l^) we obtain the result 

a- 
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In 3D the solution of the equation p^ ) at the distances r << (ap) is 

16 Q /o 1 /o\ CL 



r 



Comparing this function with the solution of the homogeneous equation (l)o{r) = C{1 — 
a/r) which is valid with the accuracy ~ a^^'^p^^'^ at r < p^^^^ we obtain the result ([TT]). 

Let us comment on the relation of our approach to that of ref.0. The function 
g{x) can be calculated using the obtained wave function of the Jastrow form. Formally 
the problem is similar to the calculation of the distribution functions of the classical 
liquid. Although the expansion in density for the corresponding partition function is 
not valid (integral / dxS{x) diverges at large distances) the formulas of ref. resulting 
in the Eq.(p^ can be obtained in the lowest order in the expansion parameter. One can 
establish the direct correspondence between the equations (|) and (|2^) . The function 
g{x) at long distances can be calculated by means of the cluster expansion (for instance, 
see 0). The connected diagrams with the lines corresponding to the function (f){x) — 1 
between the points 1 and 2 are considered. First, the Fourier transform of the function 
0(a;) — 1 may be replaced by Sk- In the integrals corresponding to the diagrams of the 
cluster expansion the long distances where the function S{x) is small are important. 
Therefore only the diagrams that does not contain the two different paths connecting 
the points 1, 2 are required. These diagrams can be refered to as a "chain" diagrams (see 
Fig.l). Adding once more line to each of these diagrams leads to an extra suppression. 
Since S{r) at r — > oo, it is clear that the sum of the "chain" diagrams gives the 
exact result for the asymptotic of the distribution function. The " chain" diagrams can 
be summed up in the momentum representation. We obtain the following relation: 

h = -r^- (24) 

Substituting the solutions of the equations (H) and (0) into Eq.(p^) we find that this 
relation is valid with the accuracy required for the estimates (0) and (|17|) . 

The same approximation can be used to calculate the Feynman structure factor 
F{k) which determines the energy of the low- lying excitations (phonons). In this case 
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the distribution function related to the square of the wave function $ (defined as in 
Eq. ([l8|) ) should be calculated. We find 

1 - 2pSk 

In contrast to the other wave functions of the Jastrow form used in the variational 
studies of the bose liquid |^ we obtain the correct behavior F{k) k at k ^ due 
to the behavior S{r) ~ l/r*^^^ at the asymptotically large distances (for a related 
discussion see ref. P]). The phonon energy 
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is in agreement with the predictions of the other approaches §]. The momentum 
distribution Uk can be readily calculated. In terms of the wave function it has the form 



nk 



N j dxe"^"^ J dx2 ■ ■ ■ dx^^ix, X2 ■ ■ ■ a:Ar)$(0, X2 ■ ■ ■ x^) 



Uk'^k = P, the normalization < $|$ >= 1 is implied). The result is as follows 

^ l-2pSk 2 \^k^ + Aapk' ) 

which coincides with the known expression 0. In two dimensions the expression 
( p!?!) should be substituted for a. In 2D the condensate is Nq = N{1 — ^a). 

Finally, the situation is different for one-dimensional problem. From the point of 
view of our method the low density limit coincide with the weak coupling limit. In this 
case the expression (^ is correct in the weak interaction limit regardless of the shape of 
the potential while for the strong coupling the perturbation theory is inapplicable which 
is in agreement with the exact solution of the problem for the 6- function potential 0. 
The same can be true for the lattice system. For instance for a system described by 



the Hamiltonian (for example, see ref. ||10||) 



^ = E (btbj + h.c.) + Uj2n,in, - 1), = bfb,, 

<ij> i 
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at U ^ and the density of order of unity Eq.(|^) is asymptotically exact in any 
dimensions. 

In conclusion, we presented the method to describe the interacting bose gas at 
zero temperature. The expansion in the irreducible functions for a logarithm of the 
ground state wave function was used. For a low density the equation for a Jastrow 
wave function was solved. It was argued that the contribution of the three-particle 
component of the wave function to the equation is suppressed. For three-dimensional 
system the leading order correction for the ground state energy in particle density was 
reproduced. For the two-dimensional dilute bose gas the ground state energy in the 
leading order in the parameter | lna^p|~^ where a is two-dimensional scattering length 
was obtained. 
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Appendix 

Here we estimate the contribution of the three-particle function S'3(xi, X2, xs) to 
the equation (j^). The three- particle equation is obtained after the integration of the 
basic equation (P over the coordinates 4, . . . A^. It is necessary to take into account the 
equations (|^) and for the terms that depend only on one of the variables xu, X13, 0:23. 
In the coordinate space the equation is 

hdl + dj + dl)S3il, 2, 3) + 9i^(12)9i^(13) + 92^(21)^2^(23) (25) 

+9i^(12)9i5(13) + . . . = (F(12) + F(13) + F(23)) 

where the dots stands for the terms of order ~ 5*25*3, 5*|. These terms are not important 
for the solution of Eq. ( pS] ) in the leading order. The right-hand side is equal to zero in 
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the infinite volume limit. In the momentum representation the function S{ki, k2, k^) is 
defined as 

J (1X1(1X2(1X36'''''''+'''^''^+'''^''^ S{1, 2, 3) = {27fy6\ki + k2 + k3)S{ki, k2, ks). 

The solution of the equation is 

2 

S{ki, k2, k-i) = - 2 I ,2 I 7.2 ihk2SkiSk2 + kiksSk^Sks + k2k3Sk2Sk3) ■ (26) 

The contribution of the function S3 to the left-hand side of Eq.(^) has the following 
form: 

2p J (lx3S{13)idid3-dl)S3{l,2,3), (27) 

which is the function of the variable 0:12 = x. Substituting the solution (|26|) to the 
integral in the momentum space corresponding to the Fourier transform of the function 
( pTl) we obtain the following expression: 

2p / , 9 (^(^ - p)SkSk^p + pip - k)SpSk-p + kpSkSp) , (28) 

Jp p^ + k^ — kp 

where k is the external momentum. It is possible to estimate the integral of Eq.(|28|) 
at different values of the momentum k with the help of Eq.([l^. For example at = 
we have 

4p / p'S'^ 
Jp 

This integral is easily estimated as O(a^) and 0{a{a^^'^p^^'^)) respectively in two- and 
three- dimensional space. Analysing Eq. (p8|) it is easy to show that these estimates are 
valid at arbitrary k. Therefore the considered contribution does not change result in 
the approximation required for the derivation of the formulas (p!TD and ([T7|). 
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